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Abstract 

We report the measurement of the charmed baryonic decay — > A^p with a branching fraction 
of (2.19t[3:^^ ± 0.32 ± 0.57) x 10"^ and a statistical significance of 5.8 a. The errors arc statistical, 
systematic, and the error of the A+ —>■ pK~iT^ decay branching fraction. This is the first observa- 
tion of a two-body baryonic B decay. The analysis is based on 78.2 fb"~^ of data accumulated at 
the T(4S') resonance with the Belle detector at the KEKB asymmetric e^e" collider. 

PACS numbers: 13.25.Hw, 14.20.Lq 
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Although the four- and three-body baryonic B decays A+j^tt+tt^ and B^ A+pvr" 

are experimentally well established 0,0]; there has been, until now, no reported observation 
of any two-body mode, such as B^ — > A'^p. In a previous Belle analysis, based on a 29.1 fb~^ 
data sample 0], we obtained the following branching fractions for four-, three- and two-body 
decays: 

^ A+:p7r+7r-) = {llMtlfj ± 1-98 ± 2.87) x 10-^ 
B{B- A+pTT-) = {l.SYt'^oil ± 0.28 ± 0.49) x 10-^ 
^ A+p) < 0.31 X 10-^ (90% confidence level). 

The measured branching fractions decrease rapidly with decreasing decay multiplicity. This 
suppression of lower multiplicity decays is a key issue in the understanding of the mechanism 
behind charmed baryonic B decays. 

There are several different theoretical calculations for the B^ A^p branching fraction, 
based on a diquark model 0, a QCD sum rule model and pole models ^,0]. They 
differ by an order of magnitude. Thus, a measurement of the two-body decay B'^ A'^p 
branching fraction would distinguish between these different theoretical approaches and 
provide important insight into the underlying physics. 

In this paper we report the first observation of the two-body decay B'^ A+p. The 
analysis is based on a data sample of 78.2 fb~^ accumulated at the T(4S') resonance with 
the Belle detector at the KEKB 8GeVe~ on 3.5GeVe"*' asymmetric collider. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50-layer cylindrical drift chamber (CDC), a mosaic of 
aerogel threshold Cerenkov counters (ACC), a barrel-like array of time-of- flight scintillation 
counters (TOF), and an array of CsI(Tl) crystals (ECL) located inside a superconducting 
solenoidal coil that provides a 1.5 T magnetic field. An iron flux return located outside the 
coil is instrumented to detect muons and i^r mesons (KLM). The detector is described in 

n 

detail elsewhere I7|l. We use a GEANT based Monte Carlo (MC) simulation to model the 

n 

response of the detector and determine its acceptance 

We detect the A+ via the A+ pK~7T^ decay channel. Inclusion of charge conjugate 
states is implicit unless otherwise stated. The particle identification (PID) information from 
the CDC, ACC and TOF is used to construct likelihood functions Lp, and L^^ for the 
proton, kaon and pion assignment for all charged tracks, respectively. Likelihood ratios 
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LR{A/B) = La/ (La + Lb) are required to be greater than 0.6 to identify a particle from A+ 
decay as type A, where B denotes the other two possible assignments among kaon, pion or 
proton. In order to maintain high efficiency for the high momentum prompt antiproton that 
comes directly from the primary meson decay, we rely more heavily on the kinematic 
reconstruction and loosen the PID requirement to LR{A/B) > 0.2, which improves the 
efficiency by a factor of about 20%. Electron and muon candidates are removed if their 
combined likelihood ratios from the ECL, CDC and KLM information are greater than 0.95. 
A candidate is selected if the invariant mass M{pK^TT^) is within 0.010 GeV/c^ (2.5 cr) 
of the 2.285 GeV/c^ A+ mass. A A+ mass constrained fit is carried out at the reconstructed 
A+ decay vertex to remove background including secondary particles from A or Ks decay. 

The B^ A+p events are identified by their energy difference AE = (J2 Ei) — -Ebeam, 
and the beam-energy constrained mass Mbc = \l E^^^^ — {J^PiY, where -Ebeam is the beam 
energy, and Pi and Ei are the three-momenta and energies of the B meson decay products, 
all defined in the center-of-mass system of the e^e~ collision. We select events with Mbc > 
5.20GeV/c^ and \ AE\ < 0.20 GeV. The prompt antiproton track and the virtual A^ track 
are required to form a common B mass constrained decay vertex. To suppress continuum 
background, we impose requirements on event-shape variables. We require | cos^thrl < 0.80, 
where ^thr is the angle between the thrust axis of the B candidate tracks and that of the 
other tracks. This requirement eliminates 80% of the continuum background and retains 
80% of the signal events. We also require R2 < 0.35, where R2 is the ratio of the second 
to the zeroth Fox- Wolfram moments jSl]. This requirement rejects 50% of the remaining 
continuum background and retains 90% of the signal. If there are multiple candidates in an 
event, the candidate with the best x\ for the B vertex fit is selected. 

Figured shows a scatter plot of AE versus Mbc and their projections for selected events. 
The AE projection is shown for Mbc > 5.270 GeV/c^ and the Mbc projection is shown for 
\AE\ < 0.030 GeV. The widths determined from single Gaussian fits to signal MC events 
are 2.7 MeV/c^ and 10.3 MeV for Mbc and AE, respectively. A two-dimensional binned 
maximum likelihood fit is performed to determine the signal yield. For this fit, the AE 
distribution is represented by a double Gaussian for the signal plus a first order polynomial 
for the background, and the Mbc distribution is represented by a single Gaussian for the 
signal plus the ARGUS function jlo| for the background. The region AE < —0.1 GeV is 
excluded from the fit to avoid feed down from modes including extra pions, which produces 
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FIG. 1: Candidate B'^ K'^p events: (a) scatter plot of IS.E versus Mbc (b) A£^ distribution 
for Afbc > 5.270 GeV/c^ and (c) Mbc distribution for [A£;| < 0.030 GeV. The curves indicate the 
result of a two-dimensional fit. 

the bump structure observed in the region Ai? < —0.15 GeV. In the fit, the signal shape 
parameters are fixed to the values fitted to the signal MC, and the signal yield and the 
background parameters are allowed to fioat. The curves in Figure ^b) and (c) indicate the 
results of this two-dimensional fit. 

The signal peak positions determined from fits to the data, (5279.5 ±0.3) M eV/ c^ for Mbc 
and (0.9 ± 1.8) MeV for Ai?, are consistent with the world average mass |ll| and zero, 
respectively. When we use single Gaussians for Mbc and /S.E signal functions and fit with the 
widths as free parameters, the fitted values in the data are found to be (1.3 ± 0.3) MeV/c^ 
and (6.9 ± 1.5) MeV, respectively, which are narrower than those determined with the signal 
MC. The probability of obtaining such narrow widths is (9(1%) and is attributed to a 
statistical fluctuation. We also investigate the decays —>■ A'^pn^TT^, B^ A^pn^ and 
B'^ ^ J / ip K* (892)'^ , J/ip pp as control samples, and find that for these modes Mbc and 
widths are consistent between the data and signal MC. The effect of the narrow widths 
to the signal yield is investigated by applying fits where single Gaussians with widths allowed 
to fioat are used for the signal shapes. The difference in the fitted yields is taken into account 
in a systematic error as discussed below. 

From the fit we obtain 19.6l|;4 signal events. From separate fits to the charge conjugate 
modes we obtain signal yields of 6.4II4 and 13.3l'l3;5 events for B^ — > A+p and B^ — >• A^p, 
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respectively. These are consistent within statistical errors. 

The branching fraction is calculated as Ns/{e x A''^^ x i3(A+ pK^n^)), using the 
measured signal yield A^s and the decay branching fraction i3(A+ —>■ pK^n^) = (5.0 ± 
1.3)% ^|. The detection efficiency e is evaluated to be 21.1% from the signal MC. The 
number of BB pairs A^^^ is (85.0 ± 0.5) x 10^. The fractions of charged and neutral B 
mesons are assumed to be the same. 

We obtain a branching fraction of 



B{BO A+p) = (2.19+°;^^ ± 0.32 ± 0.57) x 10-^ 

where the first and the second errors are statistical and systematic, respectively. The last 
error of 26% is due to uncertainty in the branching fraction i3(A^ — > pK~TT'^). 

The total systematic error of 14.8% is determined as follows. The tracking systematic 
error is estimated to be 8% in total, assuming a correlated systematic error of 2% per charged 
track, based on tracking efficiency studies with 77 — 77 and rj tt+tt^tt'^ samples. The PID 
systematic error is 10% in total, assuming a correlated systematic error of 3% per proton 
and 2% per pion or kaon, based on studies with a A ^ pn^ sample for protons; and with a 
D*^ D^Ti^, K~7T^ sample for kaons and pions. The systematic error in the fitting 

procedure and signal shape is estimated to be 7.3%, which is half of the maximum deviation 
in the branching fractions obtained with various modifications to the fitting functions: with 
a single or a double Gaussian for the AE signal, with the widths and means for both Mbc 
and AE signals fixed to MC determined values or fitted in the data. Finally, the systematic 
error in the detection efficiency due to MC statistics is 1.3%. 




M(pKV) (GeV/c^) 

FIG. 2: Invariant mass M{pK^7r^) distribution for — > A^p candidates in the B signal region. 



Figure El shows the invariant mass distribution M{pK^Tt^) for B candidates in the signal 
region \AE\ < 0.030 GeV and Mbc > 5.27GeV/c^. The curve indicates a fit result with 
a Gaussian over linear background. The fitted width of (4.0 ± 1.0)MeV/c^ and the fitted 
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mean of (2286.4 ± 1.3) MeV/c^ are consistent with values obtained from fits to signal MC 
events, which are generated assuming the world average mass We obtain a yield 
of 17.5ll\'^ events, consistent with the B signal yield mentioned above. 

We consider a contribution in the A+p signal yield from other B decays, which 

gives an uniform distribution in the Ajt invariant mass. We analyze the A+ sideband 0.015 < 
\M{pK^-K^) - M^+\ < 0.050 GeV/c^, and obtain a B signal yield of 1.2111 events, which 
is consistent with expectation of 1.4 ± 0.4 events from the B^ A^p decay MC, assuming 
our observed branching fraction. From this we estimate the other B decay contribution of 
(—O.I^q'^) events in the A+ signal region, which is negligibly small. From a simultaneous 
fit of the 5° — > A+p signal yield and the other B decay contribution in the A^ signal and 
sideband regions, we obtain a statistical significance of 5.8 a. The significance is calculated 



as ^^J —2\n{CQ/Cmax), where Cmax and Co denote the maximum likelihoods with the fitted 
s,gnal yeld and w.th the yeld fixed at zero, respecfvely. ^ 

We investigate the B'^ A^p decay in MC with all known A^ decay modes 11]. The over- 
all detection efficiency for the A'^^pK^n^ final state, including intermediate resonances, is 
found to be consistent with that calculated for the non-resonant A'^^pK'~TT~^ decay alone. 
The other A+ decays show no peaking structure in the AE and Mbc distributions. 

Finally, the overall systematics are checked by analyzing the decay mode 
B^^J/^K*{%92f observed as J/ip -> pp and K*{%92f K-n+, which contains the 
same final state particles as the mode under study. A similar analysis procedure is ap- 
plied, except for the A+ vertex fit. A J/ip is tagged if the measured invariant mass M{pp) 
is within 0.019 GeV/c^ of the J/ip mass. A K*{892)^ is tagged if the measured invariant 
mass M{K~7r~^) is within 0.2GeV/c^ of the K*(892)° mass. Vertex fits are also carried out 
at the J/ip and B vertices. We obtain a signal of 23.7l|;g events and a branching frac- 
tion of B{B^^JH)K*{m2f) = {imtoM ± 0.14) X 10"^ which is consistent with previous 



measurements I 111. 



121. 



In summary, we report the measurement of the charmed baryonic B decay B'^ —>■ A^p 
with a branching fraction of (2.19lo;4g ± 0.32 ± 0.57) x 10^^ and a statistical significance 
of 5.8 cr. This is the first observation of a two-body baryonic B decay. The branching 
fraction is found to be about an order-of-magnitude smaller than that of the three-body 
decay B^ A^pn^. This suppression is a unique feature of two-body baryonic decays; 
in contrast, the two- and three-body mesonic B decays are comparable. A pole model 0] 
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predicts a value of the branching fraction of < (1.1 — 3.1) x 10 ^ for Atp, which is 

nnn 

consistent with our measurement, while the other models p|, |5|] give substantially larger 
values. Charmless baryonic two-body decays are expected to be suppressed by an additional 
factor of I Vub/Kfep (l^. The result reported here implies that their branching fractions should 
not be much above the 10"'' level, which are consistent with the present upper limits 3]- 
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